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Dynamical coherent structure (pattern) formation in the Klein-Gordon lattice excited by peri- 
odic external field near the optical resonance is studied. It is shown that besides spatial patterns 
discovered recently (V. M.Burlakov, Phys. Rev. Lett. 80, 3988 (1998)) the spatio-temporal patterns 
can be generated in the lattice in rather broad region of excitation parameters. An influence of 
cubic anharmonicity in the lattice potential and that of random noise on the pattern formation is 
examined. 



I. INTRODUCTION 



Dynamical coherent structure (pattern) formation is a new phenomenon in dynamics of anharmonic lattices. So 
far patterns have been studied basically in continuous systems (see Ref.l and references therein) and recently in 
granular materials ■ Some soliton-like localized structures have been studied in driven anharmonic lattices ■ 
Recently it has been found that optically driven Klein-Gordon (KG) lattice with quartic anharmonicity possesses 
pattern formation, i.e. formation of a combination of few standing waves, or, in the other words, lattice spatial modes 
(LSMs), vibrating with one and the same frequency u equal to that of the driving field The important factors 
for the pattern formation are: i) modulation instability (MI) of the k = LSM directly excited by external field and 
resulting in generation of other LSMs composing the pattern; ii) destructive interference of modulation instabilities of 
different LSMs resulting in the pattern stability. Due to interference of constituent LSMs with each other the pattern 
looks like a train of equal intrinsic localized modes (ILMs) extensively studied lately ]ll|-^6| . Optical creation of a 
random ensemble of ILMs have been recently proposed by Rossler and Page for model system |l7| and by Lai and 
Sievers for realistic spin wave systems . Due to equality of the constituent ILMs the pattern formation opens a 
new way of experimental identification of ILMs based on dynamical breaking of the lattice translational symmetry. 

To make further step toward experimental identification of patterns (or ILMs) it seems expedient to study the 
pattern generation and stability in optically driven Klein-Gordon (KG) lattice with both cubic and quartic anhar- 
monicity and in the presence of thermal fluctuations (random noise). In the present paper I show that the pattern 
formation conditions in the KG lattice with cubic and quartic anharmonicity can be analyzed on a bases of the lattice 
with pure quartic anharmonicity. The pattern generation and stability regions in the space of external excitation 
parameters are presented. It is shown that at certain conditions random noise results in frustration of the stationary 
spatial patterns discovered in ||To|| giving rise to a new type of patterns. The particles vibration amplitude in the new 
patterns is modulated in both space and time so one may consider them as spatio-temporal patterns. The latter can 
be generated close to the borders of the spatial patterns stability regions and probably correspond to a transition 
regime from spatial patterns to chaos. 



II. MI OF OPTICALLY EXCITED STANDING WAVE IN THE KLEIN-GORDON LATTICE 



At the initial stage of the pattern formation the perturbation LSM with k p is amplified due to MI of the excited 
k = nonlinear LSM (referred below as to excited carrier LSM, or excited LSM). Simultaneously the higher spatial 
harmonics 2k p , 3k p , ets. will be generated due to 3- wave and/or 4- wave mixing of the excited LSM and the perturbation 
LSM. Growing up of the higher harmonics, their mutual interaction and coupling to the excited LSM finally lead to 
the pattern formation. Thus, the MI of the excited LSM works as an ignition for the pattern formation. 

Specific feature of the KG lattice is that it possesses a single optical band of vibrations. The MI of free vibrating 
waves in various anharmonic lattices including the KG one was analyzed in jp9 - 25 1 ■ Essential points for analysis of MI 
in the driven lattice are those related to i) optical excitation at frequency to which is different from the eigen frequency 
of the k — LSM; ii) phcnomcnological damping which mimic all the dissipative interactions between the phonon 
band under consideration and any other excitations in realistic system. Note, that the dissipative interactions which 
involve only the vibrations from the optical band of the KG lattice are rigorously treated in the present analysis. 

Motion equation for n — th particle of unity mass and charge in the KG lattice is 
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d 2 U n /dt 2 + 7 • dU n /dt + LU 2 U n + K 2 ■ {2U n - [/„_! - U n+1 ) + K 3 U 2 + K 4 U^ = E e^ + c.c, (1) 

where 7 is phenomenological damping constant, lu 2 , K 4 are incite and K 2 is intercite force constants, Eq is the external 
field amplitude. The trial solution of Eq.(Q) for small anharmonicity can be chosen in the form 

U n (t) = - [Vco + Vci exp (iut + i<pi) + Vc% exp [i2ut + iip 2 ) + c.c] , (2) 

where Vcj and tp are real amplitude and phase angle respectively. Two types of the KG lattice are considered. The 
lattice with pure quartic anharmonicity (LI): luq = 354, K 2 = 6 • 10 4 , K3 — 0, K4 = —3.3 • 10 6 ; the lattice with 
both cubic and quartic anharmonicity (L2): c^o = 354, K 2 = 6 • 10 4 , K3 = 6 • 10 s , K4 = 5 • 10 5 . The incite potential 
curves for a particle in the LI and L2 lattices are shown in Fig. la. A particle vibration in the right potential well 
is to be considered. One can see from Fig. la that the particle remains in the minimum until its vibration amplitude 
do not exceed the value Uf h ~ 0.2. The latter value can be considered as a natural threshold for the point defect 
formation related to the out-off minimum jump of the particle. The parameters of the LI and L2 lattices were chosen 
in a way for the Vci <x Eq dependences to be nearly the same (see Fig. lb). This choice allows to understand a pure 
influence of cubic anharmonicity not related to anharmonic shift of the optical mode frequency. For lu < luq the 
Vcj oc Eq dependences show bistability and for lu > luq these dependences are monotonous. Below only the latter 
case is considered in detail. 

Obviously Vco — Vci = in the LI lattice and one may use the rotating wave approximation (RWA) for stability 
analysis of the excited LSM. According to Rcf. [jl0| the perturbation growing rate (increment) Im(f2(<?)) can be 
determined from the equation 

lu 2 - {lu + Q) 2 - i-y • {u + fi)l • \lu 2 -{lu- Q) 2 + i-y ■ (w - = P 2 , (3) 

obtained after substitution of (^) with perturbation 

6U n — - cos{qn) ■ [Vn exp(i • {lu — Vt)t) + Vp 2 exp(i • {—lu — tt)t) + c.c] , (4) 

into Eq(]l[) (here and below the wavevector of perturbation LSM is denoted by q while that of carrier LSM is denoted 
by k). Here lu 2 — lu 2 + AK 2 sin(g/2) 2 + -i^V^, (3 = jK^V^, Vpj is the complex amplitude, and q and 17 are the 
wave vector and the complex frequency shift of the perturbation wave respectively. 

To study MI of the solution (|2|) in general case the perturbation has to be chosen in the form 



6U n — \ cos{qn) ■ [Vpo exp(— iSlt) + Vp\ exp(i • {lu — + Vp 2 exp(i ■ {—lu — £l)t)- 
V P3 exp(i • {2lu - Q,)t) + V P A exp(i • (—2u - Q.)t) + c.c] . 



(5) 



The Im(f2(q)) function then was calculated numerically equating to zero the determinant of the corresponding system 
of five liner equations. Relative increment Im(f£(g)/w) for the LI and L2 lattices and that for free carrier LSM with 
k = in the L2 lattice with 7 = are plotted in Figs.2a-c respectively for lu = 1.05wo and three values of Vci- One 
can see that the instability regions (Im(f2(q)) > 0) in the q space for the excited carrier LSM in the LI lattice are 
much broader and the instability is stronger than those in the L2 lattice. Note that the instability of free carrier 
LSM is much stronger than that of excited one (compare Fig. 2b and Fig. 2c). By symbols in Fig. 2b are shown the 
In\{rt{q) / lu) curves for the system with pure quartic anharmonicity which will be used for approximation of the L2 

lattice. The effect™ ,ua,tic tonicity constant was chosen to he K><< . K, - (l + (§)" . Vg.) - ft (ft 

and also the factor j3 in (J3J) was changed in this case to j3 e ^ = j^KiV^. The numerical factors yi and X were 
chosen for the best fit of the increment curves shown by lines in Fig. 2b. One can see that the system with the effective 
quartic anharmonicity approximates rather well the original L2 lattice if the excitation field Eq is not very high. 

This approximation allows to suitably represent the excited LSM instability regions in the {E, 7) space. The MI 
regions of the excited LSM for the perturbation wave vectors q = 7r/3 and q — tt/2 are shown in Fig. 3 by solid lines. 
The curves were obtained from Eq.(||) under conditions £l(g = tt/3) = and Q(q = n/2) = for the left and right 
closed regions respectively. Additionally, the excited LSM amplitude must be restricted by a value much lower than 
the U® h . This is because the vibration amplitude of some particles in the pattern increases compared to its value in 
the excited LSM and may result in the out-of minimum jump of these particles. In the case shown in Fig. 3 the excited 
LSM amplitude was restricted by the value 0.13 • a, i.e. about 0.7 • U° h , what cuts the MI regions for high E values. 
In fact, the difference between particles vibration amplitudes in the excited LSM on one hand and in the pattern on 
the other hand depends on damping so the MI regions can be more extended to the higher E values at high 7 values. 
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III. PATTERN SOLUTIONS FOR K 3 = AND THEIR STABILITY 



The perfect pattern can be generated starting from the seeding LSM ~ cos(g • n) at t = 0. Then under action of 
external field E = Eoe ZUJt the system will pass through two stages: a) excitation of carrier LSM; b) growing up of 
the seeding LSM due to MI of the excited carrier LSM and simultaneous generation of other LSMs due to four-wave 
mixing. Obviously, if the number of generated LSMs at the initial stage of the pattern formation is fairly restricted 
(3 -7- 5) the LSMs may grow up to large enough amplitude to be as additional nonlinear waves (carrier LSMs) in the 
lattice and to stabilize the pattern. We restrict our consideration by the total number of carrier LSMs Nlsm = 3 
and Nlsm = 4 (3-LSM- and 4-LSM patterns, respectively). The former is build up of the carrier LSMs with wave 
vectors k\ = 0, k 2 — ir/2 and k% = it (lattice constant a — 1) while the latter with those k\ = 0, k 2 = 7r/3, k$ = 2tt/3 
and &4 = 7r. The MI of the excited carrier LSM must have a maximum around g max = tt/2 in case of Nlsm — 3 and 
around g max = tt/3 for Nlsm = 4. This is the condition for Eq or, in the other words, for the excited carrier LSM 
amplitude Vc- For Nlsm — 2 (g m ax = t) no solutions have been found. 

Because of the symmetry arguments the pattern may consist of standing LSMs only. The trial pattern solutions 
are 

Uf- LSM {t) = \ [V C i exp (wt + itpx) + V C2 cos(f n) exp(iwi + iip 2 )+ 
Vc3 cos(7m) exp(iujt + if 3) + c.c] , 



Ut~ LSM {t) = \ [V C i exp (iojt + iipx) + V C2 cos(| n) exp(iu}t + i<p 2 )+ ^ 
Vc3 cos(|7rn) exp(iujt + iip^) + Vc4 cos(7rn) exp(iu>t + 1^4) + c.c] , 

for 3-LSM and 4-LSM patterns respectively. No new LSMs important within RWA will appear due to four-wave 
mixing of those chosen. Here again Vcj are real amplitudes and j£j are phase angles. According to w and 7 values a 
single stable nontrivial (all Vcj ^ 0) solution $1 of the form (|6|)-(0) was found in the region between the dotted curves 
in Fig. 3. Outside this region the solution is unstable in the sense discussed below. An additional and strongly 
unstable solution $2 exists for the 3-LSM pattern at 7 < 0.1 • ujq while for 4-LSM patterns no additional solutions 
have been found. The examples of the 3- and 4-LSM patterns of the $i-type calculated after substitution of (^|) and 
(0) respectively into Eq.(Q) with the parameters corresponding to the points a and c in Fig. 3 are presented in Fig. 4. 
Note, that indeed the patterns can be regarded as the lattices of intrinsic localized vibrations of the odd parity |]ll| , [L2| . 

Linear stability analysis of the solution (||) within RWA was given in [jl0| . Similar approach to the 4-LSM pattern 
stability requires the total perturbation to contain all perturbation waves coupled to each other via four-wave mixing, 
i.e. all spatial harmonics resulting from a product of any three carrier LSMs from Eq.(^) on a perturbation wave. 
One can see that a set of waves 

SU n = \ {exp(i • (oj - Q)t) ■ [Vpi cos(qn) + V P3 cos((| - q)n) + V P5 cos((f + q)n) + 
V P7 cos((^ - q)n) + V P9 COs((^ + q)n) + V P u cos((tt - q)n)} + 
exp(? • (— uj — Cl)t) ■ [Vp2 cos(qn) + Vpa cos((| — q)n) + Vp 6 cos(( : | + q)n) 
+V PS cos((^ - q)n) + Vp W cos((^ + q)n) + V P12 cos((tt - q)n) + c.c] } 

fulfills this condition. Indeed the set of waves @ contains spatial harmonics k p = ±q + -|m (m — 0, 1, 2, 3). After 
coupling to any three carrier LSMs from (0) it results in a set of spatial harmonics with wave vectors k p — ±q+^m±^l 
(I = 0, 1, 2, 3) which obviously can be reduced to (^). System of twelve linear equations derived after substitution of 
([?]) and (||) into Eq.(Q) was solved numerically to determine f2(q). 

The perturbation growing rate Im (17) versus wavevector is shown in Fig. 5 for 4-LSM patterns generated at the 
points a and b in Fig. 3. The pattern (a) in Fig. 5 is obviously stable, i.e. lm(il(q)) < for < q < tt, while the 
pattern (b) is unstable (Im (fl(q)) > for some q values). As some individual carrier LSMs composing the pattern are 
unstable the pattern stability can come from destructive interference between Mis of the carrier LSMs. To qualitatively 
understand this phenomenon one may imagine the simplified picture: i) perturbation LSMs of different symmetry 
compose the perturbation state Eq. (|^) via four-wave mixing; ii) this state is excited by different (by symmetry) 
coherent excitations (different carrier LSMs). Hence, whether the perturbation state will be continuously populated 
depends on the amplitudes of the coherent excitations and on the phase angle(s) between them (the influence of 
phase angle between the carrier LSMs on the pattern stability is demonstrated in Fig. 5). It may happen that at 
certain combination of the excitation parameters the perturbation state won't be populaed at all. That means the 
Ml-mediated influence of a given carrier LSM on the perturbation (JsJ) can be cancelled by that of the other(s) resulting 
in the stability of all the carrier LSMs, i.e. the pattern. 
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IV. PATTERN FORMATION IN THE PRESENCE OF NOISE 



The stability regions for both 3-LSM and 4-LSM patterns corresponding to the condition Im (f2) < (damping of the 
perturbation) are shown by dotted lines in Fig. 3. These stability regions were determined numerically via excitation 
of the pattern starting from a small amplitude seeding wave (S-patterns). Since the stability regions are inside the 
MI regions of the corresponding excited carrier LSM one might expect that the stable patterns can be spontaneously 
generated there. This is not always true, however, if the pattern is generated without any seed but in the presence 
of random noise (N-pattern). When generated from a small amplitude noise the pattern can be unstable because 
of simultaneous generation of waves with wrong (w/q is not integer) wavevectors. The examples of the N-pattern 
formation are given in Fig. 6. Numerical study was performed for 60 particle chain using the standard conservative 
scheme of numerical integration of the motion equations (|l]) with cyclic boundary conditions. Random noise of the 
value between ±0. 00005a was generated at each time step and added to the current particles positions. The panels 
notation in Fig. 6 correspond to the points notation in Fig. 3. The 4-and 3-LSM patterns generated respectively at 
the points a and c in Fig. 3 are really stable and can be attributed to the true spatial patterns since they are similar 
to those calculated theoretically (compare dotted and solid lines in Fig. 4). The pattern at the point d is inside the 
3-LSM stability region and shows periodic modulation in time similar to that of the pattern at the point b outside the 
corresponding stability region. Therefore the N-patterns at the points b and d can be regarded as spatio-temporal 
patterns rather than the true spatial patterns. 

Spatio-temporal patterns probably also possess certain stability properties though the latter are not so definite as 
for the spatial patterns. Typical temporal evolution of the spatio-temporal pattern is shown in Fig. 7. According 
to Fig. 7 the pattern corresponding to the point b in Fig. 3 is formed over ~ 30 periods of vibration and remains 
rather stable though its contrast slightly decreases in time. Relatively long-time stability of the pattern outside the 
stability region of the true spatial pattern can be qualitatively understood. Growing up of a perturbation SU n (t) in 
the presence of the unstable pattern at the point b results in the shift of the quasiharmonic resonance frequency of 
the optic LSM 

Au(k = 0) 2 = 3X4 • (5U 2 n (t)) n t 

and consequently in the increase of the frequency mismatch uj — uj(k = 0). Due to this increase the pattern approaches 
stability region or even enters it and the growing conditions for the perturbation are violated. After damping of the 
perturbation the pattern returns to the initial point and the process starts again thus giving rise to slow modulation 
of the pattern amplitude, i.e. results in formation of spatio-temporal pattern. The latter can be stable if the pattern 
modulation is perfectly periodic, i.e. the modulation amplitude doesn't grow up, otherwise it is obviously unstable. 
In the long-time scale the evolution of the unstable spatio-temporal pattern due to continues energy collapse can lead 
to a situation when the vibration amplitude of a particle exceeds the threshold value U® h what means the point defect 
formation. 



V. OPTICAL INTENSITY REQUIRED FOR THE PATTERN GENERATION 

Our numerical experiments show that to reach the threshold for the pattern formation with N^sm = 4 in the LI 
lattice the particles in the excited LSM must vibrate with the amplitude A p ~ 0.05a. The electric field strength Eq 
of the optical excitation at u> = 1.05o;o — 10 2 cm" 1 can be then estimated using motion equation in the harmonic 
approximation. Suggesting a — 5A and 7 = 0.05o>o one obtains Eq — A p m p m e [(ujQ —lo 2 ) 2 + (7a;) 2 ] 1 / 2 /(e p e e ) ~ m p /e p 
[V/cm], where m p and e p are the particle mass and charge measured in the free electron units m e and e e respectively. 
Accordingly, the field strength is of the order of 1 V/cm for light particles like electrons and of the order of 10 5 V/cm 
for ionic solids. The latter value of Eq can be reached in the laser pulses. Obviously the pulse duration T must be 
enough for the pattern to be generated and detected. The generation stage lasts over 30 — 40 periods of vibration and 
10 — 20 periods are probably needed for the pattern detection what in total means T ~ 20 ps. Hence, for the pattern 
formation in an electronic system (e.g. charge-density wave conductor) the IR laser pulses of energy W ~ 10 -11 m J 
focused into ~ 0.01 cm 2 are required. The pattern formation in an ionic system needs much higher pulse energy 
W ~ 0.1 mJ unless the smaller focusing area is used. 

Thus, in real physical experiment a system will be subjected to external field and examined during a finite time. 
Therefore it seems reasonable to use some criterion of the pattern relative stability (during finite time period) rather 
than that given by linear stability analysis. The latter, moreover, seems to be too complicated for spatio-temporal 
patterns especially for a system with realistic potential. According to the above said the pattern stability regions in 
Fig. 3 must be transformed taking into account the pattern generation and probably detection conditions (see next 
Section). 
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VI. PATTERN GENERATION AND STABILITY FOR K 3 



The 3-LSM pattern solution in this case have the form 

U*- LSM (t) = \ [Vcio + Vcu exp {iut + i<p n ) + V C12 exp (i2orf + vp 12 ) + 

cos(fn) ■ (V C 20 + V C 2i exp (iut + iip 2 \) + Vc22 exp (i2ujt + ^22)) + (9) 
cos(7rn) • (V C3 o + ^31 exp (iwf + i(p 31 ) + V C3 2 exp (i2ut + i<p 32 )) + c.c] . 

In contrast to Eqs (Q) and (Q) the pattern (Q) contains static distortion and also the higher (second) temporal harmonic 
can not be neglected. Hence the linear stability analysis of the pattern in case K 3 ^ is too complicated. Note, that 
the approximation of the system with K 3 ^ by a system with effective quartic anharmonicity gives reasonably good 
results when applied to the carrier LSM stability but not to the pattern stability. According to the estimation given in 
the previous Section, one can revise the pattern stability conditions and consider the pattern as being relatively stable 
if it is stable during few tens of periods T = 2tt/w. As a new stability criterion one may consider a requirement for 
a pattern to be stable if it shows well pronounced spatial modulation (> 10 per cent) of particle vibration amplitude 
for a time interval of about 30-50 periods of vibration. The latter time for real system is determined by the pattern 
detection time and can be even shorter. 

The new regions of the pattern relative stability were determined from numerical experiment on the bases of the 
aforementioned criterion (see shaded regions in Fig. 8). The bubbles in Fig. 8 show the pattern generation (excited 
LSM instability) regions for the LI and L2 lattices (solid and dotted lines, respectively). Dashed-dotted thin lines 
denote the envelope for the corresponding bubbles and determine the MI cut-off for the excited carrier LSMs. These 
cut-off lines are determined from the system of equations 

n( q ) = 0, 

-sr = °» (10) 

'o 2 + !^4-K 2 i-^) 2 + (7^) 2 , 

where Q(q) is determined by Eq. (H) and again K4 and must be substituted with K± and /3 e ^ if K 3 ^ 0. 
Note, that the stability regions are separated from the envelope lines because in the close vicinity to the latter the 
perturbation growing rate is too low to form the pattern with well pronounced modulation (at least 10 per cent) 
during the time of numerical experiment. 

It is important to determine the pattern generation and relative stability regions in the (E, u>) space, i.e. the space 
of external excitation parameters. They are shown by shaded regions in Fig. 9 for 3- and 4-LSM patterns in the L2 
lattice with 7 = 0.025 ■ luq. The latter value seems quite reasonable for realistic systems. The stability regions were 
obtained numerically under condition that the pattern shows more than 10 percent modulation in (U%(t)}. function 
but doesn't result in the point defect formation (particle jump from its initial minimum in Fig. la) during 120 periods 
of vibration. Even in the case of nearly fastest evolution of the 3-LSM pattern it can be considered as stable for ~ 10 
periods of vibration (see Fig. 10) though in the long time scale the system eventually arrives to chaotic state. Taking 
into account these short-living patterns one may expect the pattern formation phenomenon in rather broad region of 
the external excitation and system parameters. Note, that the particles vibration in the time interval 3 — 4 (90-120 
periods) in Fig. 10 is similar to that in the spatio-temporal pattern in Fig. 6c suggesting that the spatio-temporal 
patterns represent the intermediate state in transition from dynamical order to chaos. 

VII. CONCLUSIONS 

The dynamical coherent structure (pattern) formation was studied in the optically driven Klein-Gordon lattice 
with cubic and quartic anharmonicity in the presence of random noise. It is shown that in general the patterns are 
characterized by both spatial and temporal modulation of the particles vibration amplitude (so-called spatio-temporal 
patterns) and can be relatively stable in rather broad region of external excitation and the system parameters. 
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Figure captures 



Fig.l. a) Incite potential function for a particle in the Klein-Gordon LI- (1) and L2 (2) lattices (see the text). 

b) Vcj oc Eo plot (see Eq. 2 in the text). Curves 1, 2 and 3 correspond to j = 1 (first harmonic), 2 (second 
harmonic) and (static displacement or, zero harmonic) respectively. Lines correspond to potential of the LI lattice 
and symbols to that of the L2 lattice. The external field frequency is to — 0.9lj (solid lines) and w — 1.05w (dotted 
lines), 7 = 0.025wo- 

Fig. 2. Perturbation relative growing rate (relative increment) Im(f2(g)/u;) for lo = 1.05^0 and 7 = 0. 025a->o in the 
LI- (a) and L2 (b) lattices; (c) MI of the free (7 = 0, E = 0) nonlinear LSM with k = in the lattice L2. External 
field amplitude E was chosen for panels a) and b) in a way when the Vc\ value is equal to 0.14 (1), 0.1 (2) and 
0.07 (3). Curves (4)-(6) in the panel (b) were calculated from Eq. (||) for the lattice of Ll-type with effective quartic 



anharmonicity K\ St = K±- ( 1 + ( ^ J • Vg x ) - g ( fa j and /3 eff = j^K^V^ (see text). The E Q values for curves 



(4), (5) and (6) in the panel b) are equal to those corresponding to Vex values for curves (1), (2) and (3) respectively. 

Fig. 3. Carrier LSM modulation instability regions in the (E ,j) space enclosed by solid curves and corresponding 
to the pattern generation by external field with lo = I.Iljq (3-LSM patterns) and with lo = X.Qhuio (4-LSM patterns) 
in the LI lattice. The calculated pattern stability regions are enclosed by dotted curves. 

Fig. 4. The stable 3-LSM and 4-LSM patterns (the < U^(t) ^function) calculated (solid lines) and obtained by 
numerical simulation in the presence of noise (dotted lines) at the points a (a) and b (b) in Fig. 3. Thin dashed lines 
show the amplitude square of the carrier LSM prior the pattern formation has began. 

Fig. 5. Relative increment lm(fl(q) / lo) calculated for stable (1) and unstable (2) 4-LSM patterns (]?]) corresponding 
to the points a and b in Fig. 3. Dotted line shows the Im(f2(g)/u;) function calculated with the same parameters as 
(1) but with <pi substituted with <pi + 7r/4. 

Fig.6a-d. Molecular dynamic simulation of the N-pattern (see text) formation in 60-particle LI lattice with cyclic 
boundary conditions for 7 and Eq corresponding to the points a-d respectively in Fig. 3. Particle displacements from 
corresponding equilibrium positions are multiplied by factor 3 for (a) and (b) and by factor 2 for (c) and (d). 

Fig. 7. The < J7„(t) >t pattern obtained from numerical experiment with LI lattice at the point b in Fig. 3 and 
taken after 25 (1), 50 (2), 75 (3), and 100 (4) periods of external field vibration. The time average was taken over 10 
periods. 

Fig. 8. The pattern generation and relative stability regions (shaded) for the LI lattice (light grey) and the L2 
lattice (dark grey) for lo = l.OSwo- The bubbles enclose the 5-, 4- and 3-LSM pattern generation regions in the Ll- 
(light grey) and L2 (dark grey) lattices. Dashed lines denote the envelope for the corresponding bubbles. 

Fig. 9. The pattern generation and relative stability regions (shaded) for the L2 lattice with K2 = 30000, 7 = 
0.025luq. Solid lines enclose the pattern generation regions determined via numerical simulation for the L2 lattice 
while dotted lines enclose those obtained from equations Sl(g = tt/3) =0 and Sl(g = tt/2) = for the LI lattice with 
effective K4 and (3 (see capture to Fig. 2). In the calculations the threshold value for the carrier LSM amplitude was 
taken equal to ~ 0.1 • a. Solid symbol denotes the point at which the pattern formation process is presented in Fig. 10. 

Fig. 10. a) Molecular dynamic simulation of the N-pattern formation in the 60-particle L2 lattice for external 
field parameters corresponding to the point shown in Fig. 9. Particle displacements from corresponding equilibrium 
positions are multiplied by factor 3. Vertical dotted lines denote the starting time moments of the < U%(t) >t 
calculation shown in b). The averaging was taken over 10 periods. 
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Fig. 3. Victor M. Burlakov, Spatial- and spatio-temporal 
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Fig. 5. Victor M. Burlakov, Spatial- and spatio-temporal 
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Fig. 6. Victor M. Burlakov, Spatial- and spatio-temporal 
pattern formation in optically driven discrete systems 
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Fig. 7. Victor M. Burlakov, Spatial- and spatio-temporal 
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